Electronic transport through a double quantum dot in the spin blockade regime: 

Theoretical models 
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We analyzed the electronic transport through a double quantum dot in the spin blockade regime. 
Experiments of current rectification by Pauli exclusion principle in double quantum dots were dis- 
cussed. The electron and nuclei spin dynamics and their interplay due to the Hyperfine interaction 
were self-consistently analyzed within the framework of rate equations. Our results show that the 
current leakage experimentally observed in the spin-blockade region, is due to spin-flip processes 
induced by Hyperfine interaction through Overhauser effect. We show as well how a magnetic field 
applied parallel to the current allows excited states to participate in the electronic current and re- 
moves spin blockade. Our model includes also a self-consistent description of inelastic transitions 
where the energy is exchanged through interactions with acoustic phonons in the environment. It 
accounts for spontaneous emission of phonons which results in additional features in the current 
characteristics. We develop a microscopical model to treat the Hyperfine interaction in each dot. 
Using this model we study the dynamical nuclear polarization as a function of the applied voltage. 
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I. INTRODUCTION 

Recent transport experiments in vertical double quan- 
tum dots (DQD's) show that Pauli exclusion principle is 
important[l| in current rectification. In particular, spin 
blockade (SB) is observed at certain regions of dc volt- 
ages. The interplay between Coulomb and SB can be 
used to block the current in one direction of bias while 
allowing it to flow in the opposite one. Johnson et al.Q 
observed as well SB and spin rectification in a lateral 
DQD. Then DQD's could behave as externally control- 
lable spin-Coulomb rectifiers with potential application 
in spintronics as spin memories and transistors. 

Spin de-coherence 0, Q and relaxation processes in- 
duced, for instance, by spin-orbit (SO) scattering [B| or 
Hyperfine (HF) interaction 0] , have shown to reduce SB 
producing a leakage current. In this paper we theoret- 
ically analyze recent experiments of transport through 
two weakly coupled QD's [l[ . We consider simultaneously 
Hyperfine (HF) interaction and emission of phonons to be 
responsible of the SB lifting and the main features in the 
experimental current/ voltage (I/Vdc) curve. According 
to our calculations, inter-dot phonon-assisted tunneling 
has to be included to have a full understanding of the 
physics behind the experimental results. In the corre- 
sponding experimental, the total electron number of the 
system is fluctuating between one and two. Due to the 
different gates voltages applied between the two dots the 
left dot (n%) can have up to one electron and the right 
one {n%) can fluctuate between one and two keeping the 
sum (m + n-i) between one and two. Current flow is al- 
lowed when the electrons in each QD have antiparallel 
spins and a finite gate voltage allows one electron in the 
left dot to tunnel sequentially to the right one and fur- 



ther to the collector. However, for weakly coupled QD's 
there is a similar probability for the electron coming from 
the left lead to be parallel or antiparallel to the electron 
spin occupying the right dot. In the first case, the elec- 
tron cannot tunnel to the right dot due to Pauli exclusion 
principle and SB takes place, presenting a plateau in the 
I/Vdc curve. 



The theoretical model presented in this paper has been 
developed in the frame of rate equations. We solve self- 
consistently a system of coupled time-evolution equations 
for electronic charge occupations and nuclei polariza- 
tions. Our theoretical results reproduce the I/Vdc ob- 
served plateau due to SB and also the main current peaks. 
HF interaction is proposed as the candidate to lift SB, 
producing spin-flip (sf) of electrons and nuclei. On the 
other hand phonon-assisted tunneling is proposed, in par- 
allel with the direct tunneling, to sustain the total current 
through the device. The electrons and nuclei spin interac- 
tions brings to the Overhauser effect, which is also called 
flip-flop interaction because each time the electron flip 
the spin up to down (down to up) the nuclear spin does 
the opposite. According to measurements on QD's by 
Fujisawa et al.Q the spin-flip time, r s f > 10~ 6 s, is much 
longer than the typical tunneling time, r tun = 1 — 100 
ns, or the momentum relaxation time, r mo = 1 — 10 ns, 
meaning that spin-flip processes due to HF interaction 
are important mostly in the SB region. Our system con- 
sists of a vertical DQD under an external DC voltage in 
the presence of a magnetic field parallel to the current. 
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II. THEORETICAL MODEL 



Elastic and inelastic tunneling 



We consider a hamiltonian: 



H 



H L +H R 

H L ^y where Hl{Hr) is the hamiltonian for 
the isolated left (right) QD and is modelled as one-level 
(two-level) Anderson impurity. H^ R (HrX ) describes 
tunneling between QD's (leads and QD's)|7|) and H/ ea ds 
is the leads hamiltonian. We use a basis that contains 
20 states given by: 




|i) = |0,T);|2) = |0,|);|3) = |T,t);|4) = Ua); 
|5> = |T,i);|6) = U,t>;|7> = |o,tr>;|8> = |oar>; 

|9) = |0,U); |10) = | T,0); |11) = | 1,0); |12) = |0,T*)i 
|i3) = |0,r);|i4) = |T,r);|l5) = |l ! r); 
|i6) = |T,D;|i7) = ||,r);|i8) = |o,o) ; 

|19) = |0,|,r);|20) = |0,rn 

We have considered two levels in the right QD. Those 
states marked with (*) correspond to the excited state in 
the right QD. Double occupied states in the left QD do 
not participate in the electron transport except for high 
reverse bias. For simplicity we do not consider here this 
regime. The system of equations of motion for state oc- 
cupation probabilities p s includes scattering in and scat- 
tering out contributions fH I9L Hoj : 




FIG. 1: (Color on line). Schematic diagrams for inelastic 
phonon-assisted tunneling between weakly coupled dots, a) 
Between ground levels in each QD. b) Between the ground 
level of the left QD and the first excited level of the right QD. 
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Our neglect of coherence^ effects is appropriate for 
weakly coupled quantum dots. Wij is the transition rate 
state j to state i. 

For the contact-QD tunneling rate Wij we use a Fermi 
Golden Rule (FGR) expression, i.e., first order time- 
dependent perturbation theory 11, 12, Tj^. The expres- 
sion for the left contact-QD tunneling reads: 



W> 



L,1(1,L) 



2^ 



l + (-)larctg(^- 



7 



(2) 

where W^ i^i ^ is the tunneling probability for the left 
contact (QD) to the left QD (contact) . A similar ex- 
pression can be obtained for the the right QD to the 
right contact: 



W. 



R,2(2,K) 



-T 2 p 



1 , s 1 

- + (-)-arctg 

2 7T 



-/■tR + M2 - eVz 



2R 
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_(3) 

T is the transmission amplitude of the outer barriers, 
(T ~ 1.3 x 10- 3 meV). 7 is the width of the QD state. 
Pl,(R) i s the chemical potential in the left contact (right) 



chemical potential in the left dot (right dot). To calcu- 
late Mi, (2) we have used, apart from the corresponding 
voltage drops through the device, the intradot and in- 
terdot charging energies: 4 and 2 meV respectively [l(. p 
the two-dimensional density of states, Vl\ is the poten- 
tial drop between the emitter (left contact) and the left 
QD and V^r is the potential drop between the right dot 
and the collector (right contact). 

Inter-dot transition rates account for both elastic tun- 
neling and inelastic phonon assisted tunneling. The cor- 
responding expression for the elastic inter-dot tunneling 
is given by: 



W x . 



2(2,1) 



^1 2 2 



(Mi - M2 + eVu) 2 + 7 2 



(4) 



where T\p is the transmission through the inner barrier 
(Ti, 2 — 5x 10~ 2 meV) and V12 is the voltage drop between 
the QD's. 

For inelastic transitions, energy is exchanged with 
phonons in the environment. In other words, at T ~ 
(we have considered zero temperature in our calculations) 
the inelastic tunneling between the two dots is assisted 
by the emission of acoustic phonons, yielding a signifi- 
cant contribution to the current. This contribution has 



(ml = imeV, m_r = ^meV for zero bias[l|), Mi, (2) is the been experimentally measured by Fujisawa et al., 14\ and 
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theoretically analyzed by Brandes et al.,[l8|. In order to 
calculate the inelastic transition rate Wf 2 due to the 
emission of phonons, we have considered the theory de- 
veloped by Brandes et al., 18[. Including piezoelectric 
and deformation potential acoustic phonons the transi- 
tion rate reads: 



+ 



w d . w 
— sin — 

w w d 



(5) 



where a P i e is a piezoelectric coupling parameter 
(a pie = 0.02), e — hw = fix - u 2 + eVu, w d = c/d 
being c the sound velocity and d the distance between 
the dots. Finally 



1 



~2 



1 



U>| 7T 2 C 3 2pMC 2 h 



(6) 



where pm is the mass density and S is the deformation 
potential (3 ~ 7eV for GaAs). In Fig. 1, we represent 
schematically the inelastic contribution to / through the 
emission of phonons, between the corresponding levels of 
each QD. 



Microscopical model for Hyperfine interaction 



We calculate the electronic spin-flip scattering rate 
W^j using a microscopic model that accounts for HF 
interactions and external magnetic fields: 



H = g e p B S-B 



N, 



E 



sji + -(s+iL + s-i* + ) 



(7) 

where A is the average HF coupling constant, (A = 
90aeV for GaAs [l9|]) and / the nuclear spin . Nur) 
is the number of nuclei in the left (right) dot, (Nl = 10 
and Nr = 1.1 x 10 6 ). For simplicity we assume that 
I = 1/2. We take B to be oriented along the z direc- 
tion (current direction). The HF interaction can then be 
separated into mean-field and flip-flop contributions: 



H = H, 



H, 



where 



being, 



Hz - [geHBB + A(I z ) L{R) ]S z 



(8) 



(9) 



{Iz)l(K) = 



N, 



1 

— E (4) 



L(R) 



L(R) 



Pl(R)\Iz 



(10) 



L{R) 



N r -N l 



L(R) 



is the nuclear spin polarization 



where is the number of nuclei with spin up(down), 

in a QD. We have chosen that initially, the nuclei polar- 
ization of the left and right dots are equal to zero. 

H z has external and effective nuclear field contribu- 
tions. The latter given by: 



B,, 



On the other hand: 



H 



A 



sf 



2N 



L(R) 



MIz)l(R) 
9e^B 



^[s+p+s-ii] 



(11) 



(12) 



is the flip-flop interaction responsible for mutual elec- 
tronic and nuclear spin flips. Nuclei in similar quantum 
dots can give rise to different effective nuclear fields. This 
can be related with different effective Hyperfine interac- 
tions in each dot. An slightly distinct number of nuclei 
can explain the different Hyperfine interactions and even- 
tually the independent behavior in terms of the effective 
nuclear fields. Our model accounts for this situation with 
N L(R) and (I z )l(R)- 

Because of the mismatch between nuclear and elec- 
tronic Zeeman energies spin-flip transitions must be ac- 
companied at low temperature by phonon emission [20j]. 
Phonon absorption is not possible at low temperature. 
Therefore for B ^ Hyperfine interaction only produces 
electronic spin-flip relaxation processes. We approxi- 
mate the current-limiting spin-flip transition rate from 
parallel-spin to opposite-spin configurations by: 



1 



L(R) 



2tt 



<H sf > 



7 



(AZ e ) 



L(R) 



(13) 



where the width 7 is the electronic state life-time broad- 
ening which is of the order of peV (7 ~ 5ueV), i.e., of 
the order of the phonon scattering rate This equa- 
tion shows that a different number of nuclei or different 
splitting Zeeman can give rise to a different spin-flip rate 
in each dot. 

The splitting Zeeman is given by: 



(AZ e ) L{R) = g e u B B + —Pi 



L{R) 



(14) 



is the total electronic Zeeman splitting including the 
Overhauser shift produced by the effective nuclear B. 



(AZ 0v 



ernauser 



A 



L(R) 



-P, 



L(R) 



(15) 



We assume that a weakly coupled QD's model do no con- 
sider molecular states. The basis of states considered re- 
flects this situation. As a consequence, the exchange cou- 
pling constant is zero and the Zeeman splitting is given 
by equation (14). 
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FIG. 2: Stationary I/Vdo (B=0). At low Vdc, I takes place 
when one electron from the (1,1) spin-singlet, tunnels to the 
double occupied singlet state in the right QD (0,2). At higher 
Vdc the system reaches the states |3) = | |, j) and |4) = | j 
, J.) (inter-dot triplet states) and / drops off due to SB. The 
SB region is the plateau between the two main peaks, where 
a finite current leakage is observed due to spin-flip induced 
by Hyperfine interaction. At larger Vdc (Vdc > 6 meV) 
the chemical potential of the right lead crosses the inter-dot 
triplet state and the right QD becomes suddenly discharged 
producing a large peak in I. 



FIG. 3: (Color on line). States occupations versus Vdc for 
B = 0. In the SB region the occupation probabilities are 
dominated by states |3) = | T, T) an d |4) = | [,[)■ See blue 
(dashed) line. At larger Vdc (Vdc > 6 meV) the chemical 
potential of the drain crosses the inter-dot triplet state and 
the right QD becomes suddenly discharged. This produces 
a large peak in the current and a dramatic reduction in the 
occupation of states |3) and |4) as expected. At the same 
time other states, which are now important in the transport, 
increase their occupations. 



The expressions we propose for the the electronic spin- 
flip scattering rate W'j depend on the different pro- 
cesses: 



1,4) 
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(16) 
(17) 
(18) 
(19) 



The equations that describe the time evolution of the 
nuclei spin polarization for both dots include the flip-flop 
interaction and a phenomenological nuclear spin relax- 
ation time T re i ax w 100s [2~0 | for the scattering between 
nuclei: 



Pr. 



Pr = W 5 , 3P3 -W° f 4P4 



Trelax 

Pr 

Treiax 



(20) 
(21) 



Including spin-flip interactions, the rate equation for 
the occupation probability of |3) = | T, t) an( l |4) = | I, I) 
is: 

f>3 = W 3t ipi + W 3 jp 7 + W^iipn 

- (W 1)3 + W 7 , 3 + W n>3 + + W^) p 3 (22) 

pi = W^ 2 P2 + WifiPS + M / 04,12Pl2 

- (T4^ 2 ,4 + T^ 8 ,4 + T4^ 12:04 + T4^ 5 s/ 4 + p 3 (23) 



Total current expression 

The system of time evolution equations for the elec- 
tronic states occupations pi and nuclei polarization of 
the left and right dot is self-consistently solved. From 
that we calculate the total current through the system 
which is the physical observable of interest. The current 
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going from the left lead to the left QD is defined as: 
lL = ll+I l r = 



e{(W3,lPl + W$, 2 P2 + Wl4,12/>12 + Wl6,13Pl3) 
— (Wl,3P3 + W*2,5P5 + ^12,14^14 + W^13, l 6 Pl 6 )} T 
+e{(W 6i ipi + ^4,2^2 + Wl7,12Pl2 + Wl5,13Pl3) 

-(Wi.ePe + W2.4P4 + W 12 ,npi7 + W 13:15 p 15 )} 1 

(24) 

where the first (second) bracket, {^(O^), represents the 
up (down) current. Similar expressions can be obtained 
for the inter-dot current (/12) and for the current going 
from the right dot to the collector (Ir). In general the 
total current through the device is: 



'1,2 



Ir 



III. RESULTS 



(25) 



Experimental results []J show, for magnetic field B=0, 
a peak at low Vdc, a big plateau and a peak of high 
intensity at large Vdc, (Vdc > 6 meV). We present sim- 
ilar calculated results in Fig. 2. At low Vdc, I takes 
place when one electron from the (1,1) spin-singlet, tun- 
nels to the double occupied singlet state in the right QD 
(0,2). This process happens when the extra energy to add 
one electron to the right QD is given by a nearby gate 
voltage. At higher Vdc the system reaches the states 
3) = I T, |) and |4) = | |, |) (inter-dot triplet states) and 
I drops off due to SB. It corresponds to the observed 
plateau. This can be observed in Fig. 3, where the state 
occupation is presented versus applied bias. According 
to it, the SB region is governed by the states |3) and |4) 
with an occupation of almost 0.5 for each state. In this 
Vdc region the only way for one electron in the left QD 
to tunnel to the right one would be through an excited 
state in the right QD. Nevertheless at zero B, is too high 
in energy to participate in the transport window. De- 
spite SB, there is a finite leakage current measured in 
the SB region due to the finite probability for electrons 
in the QD's to flip their spin by interaction with nuclei. 
At larger Vdc (Vdc > 6 meV) the chemical potential of 
the right lead crosses the inter-dot triplet state and the 
right QD becomes suddenly discharged producing a large 
peak in the current. This peak is mainly due, according 
to our calculations, to the inelastic contribution to the 
current with the emission of acoustic phonons. The elas- 
tic contribution is rather small because the ground levels 
in each dot involved in the tunneling :| f , 0) — > |0, f) and 
1,0) —» 0,1), are totally out of resonance. This situ- 
ation corresponds very well with the schematic diagram 
in the Fig. la. 

A finite B parallel to the current, produces an energy 
shift experienced by the Fock-Darwin states due to its 
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FIG. 4: (Color on line). Stationary I /Vdc curve calculated 
for different B. We observe an additional peak at finite B 
which moves to lower Vdc as B increases. For different values 
of B the resonance condition: |3) = | j, T) |7) = |0, f j*} 
and |4) = | J., i) — » |8) = |0, \ [*) occurs at different val- 
ues of Vdc- A shoulder at the right side of the independent 
peak is also observed. The single line corresponds to both 
elastic transitions (direct tunneling) and inelastic transitions 
(phonon-assisted contributions). The dotted line means that 
only elastic transitions are taken into account. In this case it 
can be observed that the shoulder and the large right peak 
collapses. This demonstrates that inelastic transitions play a 
crucial role to sustain the current not only in the large right 
peak, but also in the right shoulder of the moving central 
peak. Again elastic inter-dot transitions do not contribute 
because the corresponding levels in each dot are out of reso- 
nance. The resonant peaks move to lower Vdc as B increases 
as expected. For larger B the energy shift for the excited 
state is also larger and its final position is getting closer to 
the ground state. As a result a lesser Vdc is required to get 
the resonant condition. 



coupling with the electronic orbital momentum. B cou- 
ples with the electronic orbital angular momentum of the 
first excited state of the right QD (I = 1). An additional 
shift in the electronic energy comes as well from the ad- 
ditional confinement potential induced by B [22j |. These 
B-dependent shifts and the typical energy scales of this 
problem (on-site Coulomb energy or orbital level spac- 
ing) are much larger than the Zeeman splitting. Thus, 
Zeeman splitting is neglected only in terms of tunneling 
processes. Increasing B the first excited state of the right 
QD (0, |T*) enters in the transport window and comes 
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into resonance with the ground state of the left QD. This 
opens a new transport channel, and thus I flows through 
the device and SB is lifted. In Fig. 4 we present the 
stationary I /Vdc curve calculated for different B. We 
observe an additional peak at finite B which moves to 
lower Vdc as B increases. For different values of B, the 
resonance condition occurs at different values of Vdc- I n 
Fig. 4, the single line corresponds to both elastic (direct 
tunneling) and inelastic (phonon-assisted) contributions. 
The dotted line represents only elastic inter-dot transi- 
tions. According to these results, inelastic transitions 
play a crucial role to sustain the current. They are re- 
sponsible of the large right peak and the right shoulder of 
the moving central peak. In the latter, elastic inter-dot 
transitions do not contribute because the corresponding 
levels in each dot are out of resonance (See fig. lb). The 
results presented in Fig. 4 are in good agreement with 
the experimental curve by Ono et aLy], including the 
shoulder at the right side of the independent peak. The 
oscillations however are smeared out in the experiment. 
The reason being that we did not include in our model 
damping of the bosonic system corresponding to phonon 
cavity losses |23|. It has been shown that the system is 
extremely sensitive even to very small damping |23|. 

Fig. 5 shows the charge occupation for different states 
and for the same values of B as in Fig. 4. For small B 
the electrons occupy mainly the inter-dot triplet states 
(| |, |) and | |,|)) in the SB region with a probability 
of almost 0.5 for each one. No SB removal is observed 
apart from a finite leakage current due to sf by HF inter- 
action. In this case the resonant condition between the 
ground state of the left dot and the first excited state of 
the right dot (| T,T) |0,T f) or | U) => l°> -L I*)), 
is fulfilled at Vdc > QmV. However at larger B the 
resonant condition happens inside the SB region giving 
rise to blockade lifting. The opening of this new current 
channel corresponds to a decrease in the occupation of 
3} and |4) and to an increase in the occupation of |7) 
and |8). The shoulder at the right side of the resonant 
peak produces also a removal of SB. This is due to in- 
elastic tunneling through the inner barrier assisted with 
the emission of acoustic phonons. Summarizing, SB re- 
moval is produced by three processes: the first one is 
electronic spin-flip by HF interaction. The second one is 
elastic tunneling through the right dot excited state com- 
ing into the transport window by B. Finally, additional 
inelastic contributions to the inner tunneling assisted by 
emission of acoustic phonons. 

The microscopical model for the Hyperfine interac- 
tion, allows us to study the dynamical nuclear polar- 
ization in each dot [24]]. We calculate the nuclear polar- 
izations Pl(r)i versus the applied voltage and different 
magnetic fields. We have considered that Nl — 10 6 and 
Nit = 1.1 X 10 6 . As a consequence we obtain differ- 
ent spin-flip rates for each dot. We present the obtained 
results in Figs. 6 and 7. In Fig. 6 we represent the 
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FIG. 5: (Color on line). Stationary charge occupation of the 
electronic states for the same cases of B as in Fig. 4. In the 
SB region the inter-dot triplet states(|3) and |4)) are occupied 
in the QD's, for lower values of B. For larger B the resonant 
condition between the ground state of the left QD and the 
excited state of the right QD is fulfilled inside the SB region. 
This give rise to spin blockade removal. The opening of this 
new current channel decreases pz and p4 and increases pi and 
Ps 



nuclear polarizations of left and right dot versus applied 
voltage for four different B separately B = 2,3,5 and 
8T. As we said above, for B / only electronic spin- 
flip relaxation process are possible: from spin down to 
spin up. This means that for nuclei we have the oppo- 
site process, i.e., from spin up to spin down that explains 
the negative nuclear polarization. For all magnetic fields 
studied, the right dot presents an smaller negative polar- 
ization. This is because of the larger number of nuclei in 
the right dot (see eqn. (12)). The peculiar shape that 
these graphs present can be explained if we rewrite the 
dynamical equations of the nuclear polarization of left 
and right dot. Now we have to take into account that 
only electronic spin-flip relaxation processes are allowed 
and that the nuclear spin relaxation time is very large. 
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FIG. 6: (Color on line). Nuclear polarizations of left and right 
dot versus applied voltage for four different B separately B — 
2, 3, 5 and 8T. For B ^ only electronic spin-flip relaxation 
process are possible: from spin down to spin up. For the nuclei 
we have the opposite process: from spin up to down. This 
explains the negative polarization obtained for both dots. The 
variation of Pl(R) is opposite to the electronic occupation of 
state |4) as expected. The decreasing polarization magnitude 
versus increasing B is a consequence of the splitting Zeeman 
(see eqn. 13)). 

Finally the two equations read: 

Pl = -W#p4 (26) 
PR = -W&pt (27) 

According to these expressions the calculated values 
for Pl (R) versus applied voltage are opposite to ones ob- 
tained for /?4. That means that when the occupations 
of state |4) increases (SB region) nuclear polarization 
should increase (in negative) too. However, outside of 
the SB region or when the SB is removed (new transport 
channel through state ps = |0, j J.*)) nuclear polarization 
decreases because p4 is much smaller. 

In Fig. 7, we present in two panels the nuclear polar- 
ization for both dots and all magnetic fields jointly. The 
negative values of Pl(R) decrease as the magnetic field 
increases. This happens in both dots. The explanation 



FIG. 7: (Color on line). Nuclear polarization for both dots 
and all magnetic fields jointly. The negative values of Pl(R) 
decrease as the magnetic field increases. This happens in 
both dots. If we increase B, AZ e increases and consequently 
spin-flip rate and nuclear polarization decrease as expected. 

comes readily if we observe the expression of the spin-flip 
scattering rate (eqn. 13). The highest values for this rate 
corresponds to an splitting Zeeman AZ e ~ 0, that is ob- 
tained at B — 0. In this situation spin down and spin up 
electronic states are degenerate. However if we increase 
B, AZ e increases and consequently spin-flip rate and nu- 
clear polarization decrease. Summarizing, we obtain a 
dynamical nuclear polarization that is intimately related 
with the electronic occupation through the device. 



IV. CONCLUSIONS 

In conclusion we reproduce experimental features by 
Ono et al.[l[, and show how, inelastic transitions with 
the corresponding emission of phonons play a crucial role 
sustaining the current through the device. We show as 
well that the interplay between the electron and nuclei 
spin distributions within the dots are responsible for lift- 
ing the SB yielding a leakage current. We demonstrate 
also that at finite B, the participation of excited states in 
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the current lifts SB. We develop a microscopical model to 
account for the Hyperfine interaction in each dot. Using 
this model we study the dynamical nuclear polarization 
as a function of the applied voltage. We obtain that nu- 
clear polarization is closely related with the electronic 
occupations. Our results indicate that a combination of 
B and Vdc allows to control the current through the 
device, making these systems potential components for 
spintronics and quantum computing. 
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